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After a decade of no measurements of pion and muon rare decays, PIBETA, a new 
experimental program is producing its first results. We report on a new experimental 
study of the pion beta decay, 77+ — » ir^e^u, the 7re27 radiative decay, 7r+ — > 6+^7, and 
muon radiative decay, fi — > euv^. The new results represent four- to six-fold improve- 
ments in precision over the previous measurements. Excellent agreement with Standard 
Model predictions is observed in all channels except for one kinematic region of the 77527 
radiative decay involving energetic photons and lower-energy positrons. 
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1. Motivation 

The PIBETA experimenlF at the Paul Scherrer Institute (PSI) is a comprehensive 
set of precision measurements of the rare decays of the pion and the muon. The 
goals of the experiment's first phase are: 

(i) To improve the experimental precision of the pion beta decay branching ratio, 
7r+ — > 7r°e+z^ (also referred to as TTes, or tt^), from the present ~ 4 % to ^ 0.5 %. 

(ii) To measure the branching ratio of the radiative pion decay tt — + ev^ {'^e2i, or 
RPD), enabling a precise evaluation of the pion axial- vector form factor Fa, 
and limits on the tensor form factor Ft, predicted to vanish in the Standard 
Model (SM). 

(iii) An extensive measurement of the radiative muon decay rate, fi —^ evv^, with 
broad phase space coverage, enabling a search for non- (V— A) admixtures in 
the weak Lagrangian. 

The experiment's second phase calls for a precise measurement of the tt — > ei' 
(known as 7re2) decay rate, used for normalization in the first phase. The current 
0.33% accuracy would be improved to under 0.2%, in order to provide a precise 
test of lepton universality, and, hence, of certain extensions to the Standard Model. 
In this report we focus mainly on parts (i) and (ii) above, and briefly discuss the 
first muon radiative results, (iii). 

The rare pion beta decay, 7r+ n^e'^v (branching ratio i?^^ ~ 1 x 10~*), is one 
of the most basic semileptonic electroweak processes. It is a pure vector transition 
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between two spin-zero members of an isospin triplet, and is therefore analogous to 
superallowed Fermi (S F) tr ansitions in nuclear beta decay. The conserved vector cur- 
rent (CVC) hypothesis^'"' and quark- lepton universality relate the rate of the pure 
vector pion beta decay to that of muon d ecay via the Cabibbo-Kobayashi-Maskawa 
(CKM) quark mixing matrix element Vuc^^ in a theoretically exceptionally clean 
way."^'^ Hence, pion beta decay presents an excellent means for a precise experimen- 
tal determination of the CKM matrix element Vud ] ' ' ^ ' '^ ' ^ 1 hindered only by the low 
branching ratio of the decay. 

The CKM quark mixing matrix has a special significance in modern physics as 
a cornerstone of a unified description of the weak interactions of mesons, baryons 
and nuclei. In a universe with three quark generations the 3x3 CKM matrix must 
be unitary, barring certain classes of hitherto undiscovered processes not contained 
in the Standard Model. Thus, an accurate experimental evaluation of the CKM 
matrix unitarity provides an independent check of possible deviations from the SM. 
As the best studied element of the CKM matrix, Vud plays an important role in all 
tests of its unitarity. However, evaluations of Vud from neutron decay ha ve, for the 
most part, not been consistent with results from nuclear SF decays.E^ Clearly, a 
precise evaluation of Vud from pion beta decay, the theoretically cleanest choice, is 
of interest. 

Radiative pion decay offers unparalleled access to information on the pion's 
structure. Given the unique role of the pion as the quasi- Goldstone boson of 
the strong int eraction, the implications are far reaching. In the Standard Model 
descriptioJl^ of the 7r+ e'^vj decay, where 7 is a real or virtual photon (e+e~ 
pair), the decay amplitude Ai depends on the vector V and axial vector A weak 
hadronic currents. Both currents give rise to structure-dependent terms SDy and 
SDa associated with virtual hadronic states, while the axial-vector current alone 
causes the inner bremsstrahlung process IB from the pion and positron. The IB 
contribution to the decay probability can be calculated in a straightforward man- 
ner using QED methods. The structure-dependent amplitude is parameterized by 
the vector form factor Fy [constrained by CVC to 0.0259(5)] and the axial vector 
form factor Fa that have to be extracted from experiments. 

The ratio 7 — Fa/Fv in tt ^ ev^ decay directly determines the chiral 
perturbation theory parameter sum (/g -I- Zio), or, equivalently, ue, the pion 
polarizability-I^ESEHI These quantities are of longstanding interest since they are 
among the few unambiguous predictions of chiral symmetry and QCD at low ener- 
gies. The current status of these measurements is not satisfactory, as there is con- 
siderable scatter among the various experimental determinations of Fa/Fv , and the 
accepted Particle Data Group (PDG) average value has a 14% uncertainty!^ 

Moreover, the statistical accuracy of the present experimental data on the ra- 
diative pion decaj'^'^ ^'^ cannot rule out contributions from other al- 
lowed terms in the interaction lagrangian, namely the scalar S, pseudoscalar P, and 
tensor T admixtures-l^^ Nonzero values of any of these amplitudes would imply 
new physics outside the Standard Model. In particular, reports from the ISTRA 
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collaboratioiPEl 

have indicated a nonzero tensor term, Ft — —0.0056(17). A 
careful analysis by Herczeg of the existing beta decay data set could not rule out 
such a value of Ft, which he presumed could be due to leptoquark^--'. In contrast, 
Chizhov proposed a new intermediate chiral boson with an anomalous interaction 
with matter, in order to account for the apparent non-(V— A) behavior in RPD.2H1 
Finally, radiative muon decay (RMD) provides an even better testing ground 
for non-(V— A) interaction terms than does RPD, thanks to the absence of internal 
structure in the muon. Measuring photon-positron energy distributions in RMD 
allows an evaluation of the /i"*" decay parameter fj that is predicted to be zero in 
the V-A Standard Model. The current limit on fj is loose, 0.02 ± 0.08. It is of 
considerable interest to bring the uncertainties down to a level competitive with 
the constraints from standard muon decay. 



2. Experimental Method 

The PIBETA apparatus is a large-acceptance non-magnetic detector optimized for 
detection of photons and electrons in the energy range of 5-150 MeV with high 
efficiency, energy resolution and solid angle. The main sensitive components of the 
apparatus, shown and labeled in Fig.^ arc: 




Fig. 1. Schematic cross section of the PIBETA apparatus showing the main components: beam 
entry counters (BC, ACl, AC2), active degrader (AD), active target (AT), wire chambers (MW- 
PCs) and support, plastic veto (PV) detectors and PMTs, pure Csl calorimeter and PMTs. 
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(i) beam defining plastic scintillator detectors: BC, a thin forward beam counter, 
ACi and AC2, cylindrical active collimators, AD, an active degrader, AT, a 
9-element segmented active target that stops the tt"*" beam; 

(ii) charged particle tracking and id: MWPCi and MWPC2, cylindrical chambers, 
and PV, a 20-bar segmented thin plastic scintillator hodoscope; 

(iii) a 240-element segmented spherical pure-CsI shower calorimeter, subtending a 
solid angle of 80 % of At:. 

The entire system is enclosed in a temperature-controlled Pb house lined with cosmic 
muon veto detectors. 

To collect the available TTteta decay events we recorded all non-prompt large- 
energy (above the ji evv endpoint) electromagnetic shower pairs occurring in op- 
posite detector hemispheres (non-prompt two-arm events) . In addition, we recorded 
a large prescaled sample of non-prompt single shower (one-arm) events. Using these 
minimum-bias sets, we extract tt^ and 7re2 events, using the latter for normaliza- 
tion. In a stopped pion experiment these two channels have nearly the same detector 
acceptance, and have much of the systematics in common. 

These two event classes are part a full complement of twelve fast analog trig- 
gers comprising all relevant logic combinations of one- or two-arm, low- or high 
calorimeter threshold , prompt and delayed (with respect to tt"'' stop time), as well 
as a random and a three-arm trigger, all of which were implemented in order to 
obtain maximally comprehensive and unbiased data samples. Details of the method 
are explained in detail in Ref. |57| The list of decays measured and/or used for 
normalization in our work is given in Tabled 

Table 1 . List of decays measured in the PIBETA experi- 
ment, along with the corresponding branching ratios. 



Decay Branching ratio 





fJL^V 


1.0 








~ 2.0 X 10"" 






e^u 


~ 1.2 X 10-"' 


normalize to 






~ 5.6 X 10"* 


measure 






~ 1.0 X 10-* 


measure 




77 


~ 0.9880 


measure 




e+e~7 


~ 1.2 X 10-2 


measure 




e+e^e+e" 


~ 3.1 X 10-5 






e+e~ 


~ 6.2 X 10-8 






e^uV 


1.0 


normalize to 




e+^'F7 


~ 0.014 


measure 




e+z-'Fe+e" 


~ 3 X 10-5 


measure 
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3. Results 

The building and testing of the detector components were completed in 1998, fol- 
lowed by the assembly and commissioning of the full detector apparatus. Data ac- 
quisition with the PIBETA detector started in the second half of 1999, initially at 
a reduced pion stopping rate, as planned. The experiment ran subsequently during 
2000 and 2001 at 1 MHz 7r+ stopping rate. We review below the current status 
of the analysis of the data acquired during this running period. 



3.1. Pion beta decay, 7r+ — > Tr^e+i/ 

In spite of its low branching ratio, the pion beta decay signal, marked by a nearly 
back-to-back energetic pair of neutral showers from the 7r° 77 decay, is strong 
and without appreciable background. This is demonstrated in Fig. [3 which shows 
virtually background-free histograms of 7-7 relative timing, and of the event time 
distribution following the pion stop time, <(deg). The corresponding Csl calorimeter 
energy distributions for the irp events and for the normalizing 1:^2 events are also 
given in the same figure We note the excellent agreement between the measured 
distributions and those simulated by the Monte Carlo program GEANT3. 



„ 8000 
c 

> 6000 
uj 

° 4000 - 
f 2000 - 
=^ 
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Fig. 2. Histogram of 7-7 time differences for irp decay events (dots); curve: fit with a Gaussian 
function plus a constant (top left). Histogram of time differences between the beam pion stop and 
the vr^ decay events (dots); curve: pion lifetime (top right). Csl calorimeter energy spectra for the 
pion beta decay events (bottom left) and the normalizing decay TT2e (bottom right). 



Details of the analysis method of the decay channel can be found in Refs.1281 
and EHI Applying our method and using the PDG 2004 recommended value of 
i?-P = 1.230(4) X 10-4^ US we extract the pion beta decay branching ratio: 

RlJ = [1.036 ± 0.004 (stat) ± 0.004 (syst) ± 0.003 (7ro2)] x 10"^ , (1) 
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or, in terms of the decay rate, 

r^^P = [0.3980 ± 0.0015(stat) ± 0.0015(syst) ± 0.0013(7rc2)] s'^ . (2) 

In both expressions the first uncertainty is statistical, the second systematic, and 
the third is the 1:^2 branching ratio uncertainty. Our result represents a six-fold 
improvement in accuracy over the most precise previous measuremenlPQi. Further- 
more, our result is in excellent agreement with predictions of the SM and CVC given 
the PDG recommended value range for Vu^M^ 

Rf^ ^ (1.038 - 1.041) X 10"^ (90%C.L.) , (3) 

and represents the most accurate test of CVC and Cabibbo universality in a meson 
to date. Our result confirms the validity of the radiative corrections for the pro- 
cess at the level of 4(Toxp, since, excluding loop corrections, the SM would predict 

^„o rad. corr. ^ QQ5 _ qO?) X lO^^ at 90% C.L. 

Using our measured branching ratio R^^^ , we can calculate a new value of Vud 
from pion beta decay, V,^^^^^^^ = 0.9728(30), which is in excellent agreement 
with the PDG 2004 average, ^iJ^^'"'^^ = 0.9738(5). We will continue to improve 
the overall accuracy of the 7r/3 decay branching ratio to '^0.5% by further refining 
the experiment simulation and analysis, and by adding new data. 



3.2. Radiative pion decay, 7r+ — > e+1/7 

During our 1999-2001 run we recorded over 40,000 7re27 radiative pion decays (RPD) 
events. Besides its intrinsic interest, the 'Ke2-f process is an important physics back- 
ground to other decays under study, particularly the irp decay. The pion radiative 
decay analysis has given us the most surprising result to date, and has commanded 
significant effort on our part to resolve the issue. 

The different event triggers used in our experiment are sensitive to three distinct 
regions in the 7re27 phase space: 

• region A with e+ and 7 emitted into opposite hemispheres, each with energy 
exceeding that of the Michel edge {Em — 52 MeV), recorded in the main two-arm 
trigger, 

• region B with an energetic photon {E^ > Em), and Ee+ ^ 20 MeV, recorded in 
the one-arm trigger, and 

• region C with an energetic positron (i?e+ > Em), and E^ > 20 MeV, also recorded 
in the one-arm trigger. 

Together, the three regions overconstrain the Standard Model parameters de- 
scribing the decay, and thus allow us to examine possible new information about 
the pion's hadronic structure, or non-(V— A) interactions. The RPD data are of a 
similar quality to our tt^ event set, particularly in region A, as is readily verified 
in Fig. 121 which shows histograms of e-7 time differences and event timing with 
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Fig. 3. e-7 timing difference for tt — > eiyy decay events in regions A, B, and C (left panels, top 
to bottom, respectively). The right panels plot the 7r+ — » e~^i"y event timing relative to the it~^ 
stop gate time, i(7rG), after accidental background subtraction. Monte Carlo decay functions are 
shown as full lines; best-fit values for the pion lifetime are indicated for each region. 



respect to the pion stop gate time, t(7rG). The analysis of these data is involved; 
more details can be found in Refs. 1811 and ITI 

The dependence of the region-A experimental and theoretical branching ratio 
on the value of 7 is shown in Fig. 01 (left), indicating two solutions. The positive 7 
solution is preferred by a ratio of ^ 50:1 once data from regions B and C are 
included in the analysis (right panel). We compare the experimental and theoretical 
branching ratios for the three phase space regions in Table El We note that due to 
the large statistical and systematic uncertainties present in all older experiments, 
our values are consistent with previously published measurements. The best CVC 
fit to our data yields 

7 = 0.443 ±0.015, or i^^ = 0.0115(4) with i^yEE 0.0259. (4) 

This result represents a four-fold improvement in precision over the previous world 
average Fa — 0.0116(16) It is consistent with chiral Lagrangian calcula- 
tions l^^iiSSiSSi and will lead to a correspondingly improved precision in the order 
chiral constant ho ESEi. 

Thus, our experimental 7r+ e^v^ branching ratios and energy distributions in 
kinematic regions A and C are compatible with the [V — A) interaction. The sizable 
19 % shortfall in the measured branching ratio in region B dominates the high value 
of x^/d.o.f = 25.4 (Table O, and is disconcerting. In a fit restricted to region A 
data only, we obtain 7 = 0.480 ± 0.016; this result remains unchanged if region C 
data are added to the fit. Significantly, all previous studies except on (which. 
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Minimization! Regions A+B+C 
0.443(15) 

■■••N. ^ / 



-3 -2-10 1 

7=Fa/Fv 

Fig. 4. Left plot: 7r+ — > branching ratio values as a function of 7 = Fa/ Fy. The theoretical 

parabola follows from the V — A model) . The experimental values reflect flts to region A data only. 
Right plot: minimum values of simultaneous fits to the entire data set (regions A, B, C). 



Table 2. 


Best-fit 


7r — > 


euj branching 


ratios. 


^min 

(MeV) 


(MeV) 




^cxp 

(xl0-«) 


^tho 

(xlO-») 


50 
10 
50 


50 
50 
10 


40° 
40° 


2.71(5) 
11.6(3) 
39.1(13) 


2.583(1) 
14.34(1) 
37.83(1) 




too, found an anomaly), have analyzed only data with kinematics compatible with 
our region A. An illustration of the nature of the observed discrepancy is given in 
the region- i? plots of a conveniently defined kinematic variable A in Fig. \^ Addition 
of a small negative tensor term, Ft —0.002, improves agreement with the data. 

As of this writing the PIBETA collaboration is pursuing a dedicated run to 
determine the extent and nature of the observed discrepancy more precisely. 




Fig. 5. Left panel: measured spectrum of Ai = {2Ee /rriTr) sin^ {9e-y /2) in ■7Te2-y decay for the 
kinematic region B, with limits noted in the figure. Dashed curve: three-region global best fit 
with the pion form factor Fy = 0.0259 fixed by the CVC hypothesis, Ft = 0, and F^ free. Solid 
curve: Fy = 0.0259 and = 0.0115 from the first fit, this time with Ft released to vary freely, 
resulting in Ft = —0.0018 (3). Error bars on the points at bottom of graph reflect the expected 
uncertainties in the proposed dedicated measurement. Bottom panel: same as above, but plotting 
the variable A evaluated purely on the basis of photon and positron energies: X2 = [x + y — l)/x, 
where x = 2E-,/m,r, and y = 2Ee/m,r. The two methods agree well. 
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3.3. Radiative muon decay, /x+ — > e+i/i>7 

The 1999-2001 run produced a set of some 300,000 radiative muon decay events, 
increasing the world set by two orders of magnitude. Due to an emphasis in the early 
analysis on the pion rare decay channels, the RMD analysis has not yet reached the 
same sub-1 % level of precision. The analysis is currently at the ~ 1 % level and both 
integrated and differential branching ratios are in good agreement with the V— A 
SM predictions at this level.l^ Of particular interest are the shapes of differential 
distributions of A, a suitably defined kinematic variable (analogous to A in RPD), 
shown in Fig. |B1 Clearly, nothing like the ~ 20 % discrepancy of region B in RPD 
is observed here. 




Fig. 6. Dots: distributions of the kinematic variable A = 1 — pe cos6e-y / Ee for radiative muon 
decay events collected with the two-arm (left) and one-arm trigger (right). Histograms: Standard 
Model V-A GEANT simulation. 



This work is in progress, and the 2004 dedicated RPD run is expected to more 
than double the RMD event set, adding high-quality new data. Interpretation of 
the PIBETA RMD data below the 1 % accuracy level will not be possible without 
new reliable radiative corrections for this process. 

4. Summary 

Over a decade after the last precision measurements of the rare pion and muon 
decays were completed, the PIBETA experiment at PSI is revisiting the field, and 
producing new results. Unlike any experiment exploring the field before it, the 
PIBETA project measures simultaneously practically all of the rare tt and fi decays. 
This gives it an unusually powerful set of built-in consistency checks. 

The pion beta decay precision has been improved six-fold, resulting in a first 
definitive test of the CVC and radiative corrections in a meson. Work on this decay 
channel is continuing, and will reach 0.5 % accuracy in this first phase of the 
project. The current result, with a combined systematic and statistical uncertainty 
of 0.64% is in excellent agreement with the SM predictions. 

The PIBETA radiative pion decay measurements cover a broader region of phase 
space than previous experiments, with more than an order of magnitude higher 
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statistics, and have brought about a four-fold improvement in the precision of the 
pion axial form facto. However, the unexpected and pronounced deviation from the 
V— A description of the process is limiting this precision, and raising important 
questions. This matter is being addressed in a current dedicated run. 

Comparable improvements in precision are expected in radiative muon decay. 

Finally, in the following, second phase of the project, the PIBETA collaboration 
will turn its attention to the n ^ ev process which provides the best test of lepton 
universality and a selective check on possible new physics. 

This material is based upon work supported by the National Science Foundation 
under Grants No. 0098758 and 0354808. 
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